(N 

O 
(N 

o 

Q: 
a^: 



d 



Uncertainties in nuclear transition matrix elements for /?^/3+ and e(3^ modes of 
neutrinoless positron double-/? decay within PHFB model 

P. K. Rathi, R. Chandra^, K. Chaturvedi^, P. Lohanii, P. K. Raina"* and J. G. Hirsch^ 
^Department of Physics, University of Lucknow, Lucknow-226007, India 
^Department of Applied Physics, Babasaheb Bhimrao Ambedkar University, Lucknow-226025, India 

^Department of Physics, Bundelkhand University, Jhansi-284128, India 
'''Department of Physics, Indian Institute of Technology, Ropar, Rupnagar - I4OOOI, Punjab, India 
^Instituto de Ciencias Nucleares, Universidad Nacional Autonoma de Mexico, 04510 Mexico, D.F., Mexico 

(Dated: December 21, 2012) 

Uncertainties in the nuclear transition matrix elements M'*^"' and M'*^^' of the double-positron 
emission (/?''' /3^)oiy and electron-positron conversion (e/3^)oi/ modes due to the exchange of light 
and heavy Majorana neutrinos, respectively, are calculated for ^^Ru, ^°'^Pd, ^"^Cd, ^^'*Xe, ^"^"Ba 
and ^^^Dy isotopes by employing the PHFB model with four different parameterization of the 
pairing plus multipolar two-body interactions and three different parameterizations of the Jastrow 
short range correlations. In all cases but for ^^'^Ba, the uncertainties are smaller than 14% for light 
Majorana neutrino exchange and 35% for the exchange of a heavy Majorana neutrino. 
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I. INTRODUCTION 

The Majorana nature of the neutrinos could be imme- 
diately established by confirming the possible occurrence 
of any one out of four experimentally distinguishable 
modes of lepton number violating neutrinoless double 
beta (/3/3)o^ decay, namely the double-electron emission 
(/3~/3~)o,y, double-positron emission (/3+/?+)o,y, electron- 
positron conversion (e/3+)oj/ and double-electron capture 
(e£)oi/. The latter three modes are energetically compet- 
ing and we shall refer to them as (e"*"/3/3)p^ decay. The 
kinetic energy release in the {ee)^^ mode is the largest. 
However, the conservation of energy-momentum requires 
the emission of an additional particle in the {ee)^^ mode. 
The absorption of atomic electrons from the fC-shell is 
forbidden for the 0+ — O"*" transition due to the emission 
of one real photon. Consequently, various processes such 
as internal pair production, internal conversion, emission 
of two photons, L-capture etc. [l| have to be considered. 
The decay rates of the above mentioned processes have to 
be calculated at least by the third order perturbation the- 
ory and are suppressed by a factor of the order of 10^^ in 
comparison to the (e/3"'")oiy mode. Hence, the experimen- 
tal as well as theoretical studies of (e"'"/3/3)oi/ decay had 
been mostly restricted to (/3"'"/3"'")oi/ and (e/3"'")oi/ modes 
only. 

The idea behind the resonant enhancement of (e£)oi/ 
mode has been recently reanalyzed [1, 0| and it 

has been shown that there will be resonant enhancement 
of the (ee)o!/ mode upto a factor of 10^ provided the nu- 
clear levels in the parent and daughter nuclei are almost 
degenerate i.e. Q — {E2P — £23) ^ 1 keV, where the 
energy difference is for atomic levels. Subsequently, de- 
tailed theoretical studies on the resonant enhancement of 
(e£)o!y mode have also been performed In the mean 

time, experimental studies on resonance enhancement of 
{ee)ou mode in ^^Qe [ll|ll|, ^^Rujll, ^06q^ [ll,|i|, 
"2Sn dMi, ^^'^Ce [U and ^^ow [2lf isotopes have al- 



ready been carried out and the study of this (ee)oi/ mode 
is emerging as an interesting possibility for the investiga- 
tion of (/3/3)g^ decay. 

In addition to establishing the Dirac or Majorana na- 
ture of neutrinos, the observation of (/3/3)q^ decay can 
also ascertain the role of various mechanisms in different 
gauge theoretical models [l^]. The study of {PP)^^ decay 
can clarify a number of issues, such as the origin of the 
neutrino mass, their absolute scale as well as hierarchy, 
and possible CP violation in the leptonic sector. The 
(/3/3)qj, and {e'^ filS)^^ decay modes can provide us with 
similar but complementary information. The observation 
of (e''"/3/3)oiy decay modes would be helpful in determining 
the presence of mass mechanism or right handed currents 
[231 . The varied scope and far reaching nature of the ex- 
perimental and theoretical studies on the (/3/3)n., decay 
have been recently reviewed by Avignone et al. [24| , Ver- 
gados et al. and Faessler et al. [2^ 

The nuclear /3/3 decay proceeds through strongly sup- 
pressed channels which are very sensitive to details of the 
wave functions of the parent, intermediate and daugh- 
ter nuclei. Hence, the calculations of non-collective nu- 
clear /3/3 decay related observables are quite challeng- 
ing. In any nuclear model, there are three basic ingredi- 
ents, namely the model space, the single particle energies 
(SPEs) and the effective two body interactions. Usually, 
these are chosen on the basis of practical considerations. 
While all models are able to reproduce most of the ob- 
served {(3/3)^^ decay half lives by adjusting free parame- 
ters in the model, different predictions are obtained for 
other observables, like the (/3/3)o^ decay half lives, due to 
the inherent freedom in choosing the basic ingredients of 
the model. 

A variety of nuclear models is currently employed in 
this endeavor. Large scale shell model calculations are 
quite successful |27h29| | , but highly limited in the descrip- 
tion of medium and heavy mass nuclei. The most popular 
and successful model is the Quasiparticle Random Phase 
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Approximation (QRPA) and its extensions [30|,|3l|. The 
inclusion of nuclear deformation has also been carried out 
in the deformed QRPA [3^ l33|. the Projected Hartree- 



Fock-Bogoliubov (PHFB) '34U365, the pseudo-SU(3) .37], 
the Interacting Boson Model (IBM) ^38j , and the Energy 
Density Functional (EDF) [331 approaches. In the study 
of both {(3/3)21, ^^'^ ^Pl^)ou decay modes, the renormal- 
ized value of axial vector coupling constant qa is a major 
source of uncertainty. In the {(3/3)^^ decay the role of 
pseudoscalar and weak magnetism terms [io, El is cru- 
cial, and the finite size of nucleons (FNS) and short range 
correlations (SRC) play a decisive role vis-a-vis the radial 
evolution of nuclear transition matrix elements (NTMEs) 

[HlsllsElil. 

Usually, three different approaches have been adopted 
for estimating the uncertainties in NTMEs for (/3~/3~)gt, 
decay. The spread between all the available calculated 
NTMEs has been used as the measure of the theoreti- 
cal uncertainty [i^. The same spread between NTMEs 
can also be translated into average and standard devia- 
tion, which can be interpreted as theoretical uncertainty 
[il . [isj . According to Bilenky and Grifols [Hi , the ob- 
servation of (/5/3)gj^ decay of different nuclei will provide 
a method, in which the ratios of the NTMEs-squared can 
be compared with the ratios of observed half-lives T^/'2 
and the results of calculations of NTMEs can be checked 
in a model independent way. 

The theoretical uncertainties were estimated by Rodin 
et al. [l^l by considering two models, QRPA and 
RQRPA, with three sets of basis states and three realis- 
tic two-body effective interactions based on the charge 
dependent Bonn, Argonne and Nijmen potentials. It 
was found that the variances were substantially smaller 
than the average values and the results of QRPA, albeit 
slightly larger, are quite close to the RQRPA values. The 
critical analysis of the advantages and defixnencies in the 
approach of Rodin et al. [l^l by Suhonen 48] and Rodin 
et al. [49|] is quite instructive. Further studies on the 
uncertainties in NTMEs due to SRC using the unitary 
correlation operator method (UCOM) |53| and by self- 
consistent coupled cluster method (CCM) |5lj have also 
been carried out. 

Recently, the uncertainties in the (/3~/3~)p,, NTMEs 
due to the exchange of light 35] and heavy [3g Majorana 
neutrinos have been calculated in the PHFB model by 



employing four different parameterizations of the pairing 
plus multipolar effective two body interaction and three 
different parameterizations of Jastrow type of SRC. In 
the present work, we employ the same formalism for es- 
timating uncertainties in NTMEs for (e/3+)oi/ and (ee)oi/ 
modes of 96Ru, i02pj^ loe^^j^ i24xe, isoga and ^^^Dy 

isotopes for the 0+ — > 0+ transition. The article is or- 
ganized as follows. A brief discussion of the theoretical 
formalism is presented in Sec. II. In Sec. Ill, we ana- 
lyze the role of the different parameterizations of the two 
body interaction, the finite size of nucleons and higher 
order currents (HOC). The influence of the SRC in the 
radial evolution of the NTMEs is also presented. In the 
same Sec. Ill, we estimate the uncertainties, which are 
subsequently employed for extracting bounds on the ef- 
fective mass of light neutrinos (m^) and heavy neutrinos 
{Mn)- In Sec. IV, the conclusions are presented. 



II. THEORETICAL FORMALISM 

In the Majorana neutrino mass mechanism, the half- 
lives T^J^ for the 0+ 0+ transition of (/3+/3+)o^ and 
(e/3+)g^ modes are given by P, SO, El 



T^j2 = Goi (/3) 



(1) 



rUe (Mn) 
Here, (3 denotes the (/?'*' /?^)o,y / (^/3^)oi/ modes, 

(m^) = V' C/2 m„ m, < 10 eV, (2) 

^ — 

and 



> 1 GeV. (3) 
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where K — Ov (ON) denotes the exchange of light (heavy) 
Majorana neutrino mechanism. 

In the PHFB model, the NTMEs M(^) for the 
(/3+/3+)q^ and (e/3^)oy modes are calculated by employ- 
ing the closure approximation j34| 
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where 
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with 



(7) 



r 



and the expressions for 



n(z N) (Z-2.N+2)(S), fz,N 



3J|. The three com- 



and Fz^N{d) are given in Ref. 
ponents of the nuclear transition matrix element Af*^^^ 
are denoted by F, GT and T corresponding to Fermi, 
Gamow- Teller and tensor terms. 

The neutrino potentials due to the exchange of light 
and heavy neutrinos between nucleons having finite size 
are given by 



where fa{qrnm) = jo {qrnm) for a = F,GT and 
fa {qrnm) = 32 {qrnm) for « = T. The above expres- 
sions for the NTMEs M^^' were obtained by including 
pseudoscalar and weak magnetism terms in the nucleonic 
current and employing the Goldberger-Treiman PCAC 
relation for the induced pseudoscalar term [40] • 
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a \ nm } 



2R f fa{qr nrn ) 



TT J {q + A) 
2R 



{mpme)n 
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Usually, the influence of the finite size of nucleons 
(FNS) is taken into account through dipole form factors. 
The functions /if (9), horiq) and hT{q) are written as 
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where 

9v{q') 
9A{q') 

9M{q^) = 

9p{q^) = 

with gv — 1.0, gA 
0.850 GeV and 
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= 9A 
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1.254, K ^ ^p — fin = 3.70, Av — 
1.086 GeV. 

Consideration of Eq. ©-Eq. (HI]) and Eq. © im- 
plies that the Fermi matrix element m'^j^^ has one term 



-g\M^p^Jyy, the Gamow- Teller matrix element 

has four terms, namely M^_^^, ^'^gt~ap^ ^gt-pp^ 

riK) 
-ppi 



MM ^^'^ there are three terms M!p^\p, 
^T^MM associated with the tensor matrix element 



In the literature, the short range correlations (SRC) 
have been included through the exchange of w-meson [svj , 
effective transition operator ^55] , unitary correlation op- 
erator method (UCOM) 42, 50], self-consistent CCM [5j| 
and phenomenological Jastrow type of correlations with 
Miller-Spenser parameterization [53|. Further, Simkovic 
et al. [51] have shown that in the self-consistent CCM, 
it is possible to parametrize the effects of Argonne V18 
and CD-Bonn nucleon-nucleon (NN) potentials by the 
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Jastrow correlations with Miller-Spenser type of param- 
eterization given by 

/(r) = l-ce-"'''(l-6r2). (14) 

In the present work, the above form is adopted with a = 
1.1 /m-2, 1.59 1.52 /m~^ b = 0.68 /m^^ 1.45 

/m-^ 1.88 and c = 1.0, 0.92, 0.46 for Miller- 

Spencer parameterization, Argonne V18 and CD-Bonn 
NN Potentials, which are denoted as SRCl, SRC2 and 
SRC3, respectively. 

The NTMEs M'-^' of the (/3+;3+)o^ / (e;S+)o^ decay 
mode in the PHFB model have been already discussed 
in Ref. js^. The same formalism is employed here. The 
axially symmetric HFB intrinsic state |<i>o) with K = 
specified completely by the amplitudes {uim,Vim) and 
expansion coefficients Cij^rm is obtained by minimizing 
the expectation value of the effective Hamiltonian given 
by H 

H = Hsp + V{P) + ViQQ) + V{HH), (15) 

in a basis constructed by using a set of deformed states. 
Here, Hgp denotes the single particle Hamiltonian and 
V{P), V{QQ) and V{HH) are the pairing, quadrupole- 
quadrupole and hexadecapole-hexadecapole parts of the 
effective two-body interaction, respectively. 

The details about the parameters of the pairing force 
Gpp and G„„ as well as three strength parameters of 
quadrupolar interaction, namely the proton-proton X2pp, 
the neutron-neutron y-inn and the proton-neutron X2pn 
have been given in Refs. [33,[5l,[5l]. Specifically, X2pp — 
X2nn = 0.0105 MeV6~^, where h is the oscillator parame- 
ter and the strength parameter X2pn was varied to fit the 
experimental excitation energy of the 2+ state, i?2+- 
Presently, we employ in addition an alternative isoscalar 
parameterization by taking X2pp = X2nn = X2pn/2 and 
the three parameters are varied together to fit E2+ ■ 
These two parameterizations of the quadrupolar interac- 
tion are referred as PQQl and PQQ2. The details about 
the HH part of the effective interaction V{HH) have 
also been given in Ref. [l^. The calculations including 
the hexadecapolar term HH are denoted as PQQHH . 
With the consideration of the hexadecapolar interaction, 
we end up with four different parameterizations, namely 
PQQl, PQQHHl, PQQ2 and PQQHH2 of the effec- 
tive two-body interaction. By employing the four differ- 
ent parameterization of the two body effective interac- 
tion and three different parameterizations of SRC, sets 
of twelve NTMEs Af^O"^) and Af(o^) for the (/3+/3+)o^ 
and (£/3^)oiy rnodes are obtained using Eq. ([SJ and sub- 
sequently, the mean and standard deviations are calcu- 
lated for estimating uncertainties associated in the results 
of the present work. 

III. RESULTS AND DISCUSSIONS 

In the present work, we use the same model space, 
single particle energies (SPE's) and parameters of the ef- 



fective two-body interaction as our earlier calculations 
on the (e+/3/3)2^ decay of ^^Ru, i02pd, W6,W8(j^ 
I24,i26xe, I30,i32gj^ jggj I56j3y [5^ isotopcs for the 
0+ — > 0+ transition. The calculated yrast spectra, re- 
duced B{E2:0'^ — ^ 2+) transition probabilities, static 
quadrupole moments Q(2+) and gyromagnetic factors 
g{2~^) j55l457| are in an overall agreement with the exper- 
imental data due to PQQl parameterization. The maxi- 
mum change in all the calculated spectroscopic properties 
is about 10% except for the B{E2:0+ -J> 2+) and g{2+) 
of ^-^^Ba, which change by 19.5% and 16.6%, respectively, 
employing the other three parameterizations. 

A. Effects due to finite size of nucleons and sliort 
range correlations 

The theoretically calculated sets of twelve NTMEs 
M^""^) and M'"^^ using the HFB wave functions in con- 
junction with PQQl, PQQHHl, PQQ2 and PQQHH2 
interaction and three different parameterizations of the 
Jastrow type of SRC for ^^Ru, ^o^Pd, ^°^Cd, ^"^^Xe, "OBa 
and ^^^Dy nuclei are given in Table ID The sets of twelve 
NTMEs Af(°'') and M(°^) are calculated in the approx- 
imation of finite size of nucleons with dipole form factor 
(F) and finite size plus SRC (F+S). Further, the NTMEs 
^(Of) g^j-g calculated for A and ^/2 in the energy denom- 
inator in the case of F-l-S. We present the components of 
NTMEs M(°'') as weU as M'-°^\ namely Fermi, Gamow- 
Teller and tensor matrix elements of ^"^Cd in Table |ll] 
for explicitly displaying the role of higher order currents 
(HOC). In Table ini the changes in NTMEs M^^"^ and 
]Vf(OJV)^ due to different approximations are displayed. 
The following observations are noteworthy. 

(i) Changing A to A/2 in the energy denominator, the 
changes in the NTMEs M^O") vary between 9%- 
12 % exhibiting that the dependence of NTMEs on 
average excitation energy A is small, which supports 
the use of the closure approximation in the case of 
the (/?/?) 0^ decay. 

(ii) Inclusion of effects due to FNS induces changes in 
the NTMEs ^^^"^^-hM^"^^ by 9.0%-13.0%. Further, 
the addition of higher order currents (HOC) reduces 
the NTMEs by 11.0%-15.0%. 

(iii) With the addition of SRCl, SRC2 and SRC3, the 
NTMEs M'-^"'^ vary approximately by 13%-20%, 
0.9%-2.6% and 2.5%-3.0%, respectively, in compar- 
ison to the case F. 

(iv) The NTMEs M^"^^ + M^^P in the case of heavy 
neutrino exchange, with the consideration of FNS 
instead of point nucleons, vary by 18.0%-29.0%, and 
the inclusion of HOC results in further reduction by 
about 14.0%-17.5%. 

(v) In the case of heavy neutrino exchange, the NTMEs 
]{f{0N) bgcome smaller by approximately 65%- 
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TABLE I: Calculated NTMEs M^O"' and M'"^' in the PHFB model with four different parameterization of effective two-body 
interaction, namely (a) PQQl, (b) PQQH HI, (c) PQQ2 and (d) PQQH H2 and three different parameterizations of Jastrow 
type of SRC for the {P+P^)^^^ and {eP+)^^ modes of ^^Ru, i°2pd, ""^Cd, ^^^Xe, "'^Ba and ^^^^Dy isotopes due to the exchange 
of light as well as heavy Majorana neutrinos. See footnote on p. 3 of Ref ^35*] for further details. 

Nuclei Light neutrino exchange Heavy neutrino exchange 



F F+S F+S(A/2) F F+S 











SRCl 


SRC2 


SRC3 


SRCl 


SRC2 


SRC3 






SRCl 


SRC2 


SRC3 


«'=Ru 


(a) 


4 


7979 


4 


1538 


4 


7331 


4 


9191 


4 


5739 


5 


1844 


5 


3777 


251 


3551 


87 


5090 


150 


6483 


205 


2102 




(b) 


4 


7820 


4 


1352 


4 


7164 


4 


9032 


4 


5501 


5 


1626 


5 


3567 


252 


6125 


88 


1241 


151 


5084 


206 


2858 




(c) 


4 


8334 


4 


1861 


4 


7686 


4 


9555 


4 


6107 


5 


2245 


5 


4187 


252 


8381 


88 


1649 


151 


6459 


206 


4813 




(d) 


4 


7399 


4 


1000 


4 


6753 


4 


8601 


4 


5120 


5 


1182 


5 


3103 


250 


2061 


87 


4415 


150 


1862 


204 


3879 


102p^ 


(a) 


5 


3695 


4 


5877 


5 


2981 


5 


5230 


5 


0512 


5 


8007 


6 


0346 


296 


4236 


97 


6551 


174 


1164 


240 


2760 




(b) 


4 


5052 


3 


8203 


4 


4407 


4 


6377 


4 


1853 


4 


8400 


5 


0449 


261 


2562 


87 


2707 


154 


2305 


212 


1479 




(c) 


5 


4006 


4 


6167 


5 


3292 


5 


5547 


5 


0831 


5 


8348 


6 


0694 


297 


6968 


98 


3279 


175 


0482 


241 


4044 




(d) 


4 


4595 


3 


7823 


4 


3959 


4 


5907 


4 


1432 


4 


7907 


4 


9933 


258 


6259 


86 


5570 


152 


7985 


210 


0766 


i"«Cd 


(a) 


8 


4560 


7 


2607 


8 


3403 


8 


6835 


8 


0547 


9 


1947 


9 


5519 


452 


6855 


149 


5474 


265 


9923 


366 


9093 




(b) 


6 


9410 


5 


9037 


6 


8370 


7 


1347 


6 


5165 


7 


5021 


7 


8119 


394 


7635 


132 


0821 


233 


0127 


320 


4765 




(c) 


8 


5399 


7 


3370 


8 


4229 


8 


7683 


8 


1436 


9 


2902 


9 


6497 


455 


6637 


150 


6151 


267 


7800 


369 


3362 




(d) 


7 


7425 


6 


6175 


7 


6293 


7 


9524 


7 


3228 


8 


3911 


8 


7274 


428 


1103 


143 


0153 


252 


5179 


347 


4477 


Ae 


(a) 


4 


1442 


3 


5405 


4 


0770 


4 


2507 


3 


9471 


4 


5153 


4 


6966 


230 


OOlO 


76 


6774 


135 


6415 


186 


8959 




(b) 


3 


4015 


2 


8367 


3 


3342 


3 


4963 


3 


1370 


3 


6639 


3 


8331 


213 


1788 


70 


0146 


124 


7059 


172 


4178 




(c) 


3 


6899 


3 


1428 


3 


6275 


3 


7849 


3 


5024 


4 


0157 


4 


1799 


207 


7699 


68 


5459 


121 


8099 


168 


1977 




(d) 


3 


4722 


2 


8994 


3 


4045 


3 


5690 


3 


2056 


3 


7406 


3 


9123 


216 


5439 


71 


2281 


126 


7645 


175 


1900 


""Ba 


(a) 


3 


5986 


3 


0605 


3 


5369 


3 


6914 


3 


4204 


3 


9254 


4 


0868 


205 


5885 


68 


6378 


121 


1524 


166 


7736 




(b) 


2 


8901 


2 


4039 


2 


8372 


2 


9769 


2 


6473 


3 


1067 


3 


2525 


186 


0534 


62 


0004 


109 


6716 


150 


9816 




(c) 


2 


9496 


2 


4910 


2 


8950 


3 


0266 


2 


7799 


3 


2083 


3 


3457 


174 


2028 


57 


7641 


102 


3282 


141 


1255 




(d) 


1 


5194 


1 


2183 


1 


4801 


1 


5662 


1 


3358 


1 


6133 


1 


7032 


110 


7898 


34 


9642 


63 


7693 


89 


0510 


l56Dy 


(a) 


2 


1901 


1 


9136 


2 


1712 


2 


2513 


2 


1304 


2 


4044 


2 


4882 


112 


6022 


39 


8413 


68 


3273 


92 


5002 




(b) 


1 


9284 


1 


6658 


1 


9097 


1 


9857 


1 


8415 


2 


1010 


2 


1805 


106 


1906 


37 


2359 


64 


1897 


87 


1017 




(c) 


2 


5208 


2 


1991 


2 


4982 


2 


5913 


2 


4520 


2 


7702 


2 


8675 


130 


9372 


46 


4436 


79 


5200 


107 


5918 




(d) 


2 


5789 


2 


2241 


2 


5526 


2 


6551 


2 


4674 


2 


8168 


2 


9241 


142 


7229 


49 


7900 


86 


0621 


116 


9468 



TABLE II: Decomposition of NTMEs M^^'^^ and M(°'^' for the {P'^P'^)^^ and (£^5+),,^ modes of ^°'^Cd including finite size 
effect (F) and SRC (F+S) for the PQQl parameterization. 



NTMEs 


Lif 


jht neutrino exchange (K 


= Qv) 


Heavy neutrino exchange [K = 


QN) 




F 




F+S 




F 




F+S 








SRCl 


SRC2 


SRC3 




SRCl 


SRC2 


SRC3 




2.6164 


2.2956 


2.6164 


2.7064 


133.9749 


69.7853 


105.3773 


125.4157 


^"GT-AA 


-7.8995 


-6.7491 


-7.7996 


-8.1279 


-461.3370 


-179.9510 


-301.7550 


-393.8500 


^"GT-AP 


-.5562 


-.3643 


-.5003 


-.5570 


218.6080 


55.4365 


118.4700 


172.6570 


^"gt-pp 
M^gV 


1.8965 
-.2945 
-6.8537 


1.3980 
-.1479 
-5.8634 


1.7872 
-.2276 
-6.7404 


1.9323 
-.2735 
-7.0261 


-83.4189 
-52.2544 
-378.4023 


-11.3250 
20.1072 
-115.7323 


-36.4994 
9.1939 
-210.5905 


-60.7425 
-16.8909 
-298.8264 




-.0306 


-.0311 


-.0319 


-.0318 


14.8474 


14.5508 


15.8270 


15.8681 


]S/Lrj-\ P P 


.0835 


.0849 


.0866 


.0864 


-5.9893 


-5.8221 


-6.4169 


-6.4410 




.0087 


.0087 


.0092 


.0091 


2.0565 


1.8336 


2.2001 


2.2449 




.0615 


.0625 


.0638 


.0637 


10.9146 


10.5624 


11.6101 


11.6720 




8.4560 


7.2607 


8.3403 


8.6835 


452.6855 


149.5474 


265.9923 


366.9093 
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TABLE III: Changes (in %) of the NTMEs M^°'''> and A'/'°^' due to exchange of hght and heavy Majorana neutrinos, respec- 
tively, for the {P'^ S'lid {sP^)^^ modes with the inclusion of finite size effect (FNS) as well as finite size effect -l-HOC (F), 
F+SRC (F+SRCl, F+SRC2 and F+SRC3) for four different parameterizations of the effective two-body interaction. 







Light neutrino exchanj 


;e 




Heavy neutrino exchanj 


;e 


FNS 


F 


F+S 


FNS 


F 


F-fS 


(a) 


9.08-10.31 


10.91-12.47 


12.63-14.95 


25.89-28.91 


14.41-17.51 


64.62-67.06 








0.86-1.72 






39.32-41.26 








2.53-2.86 






17.85-18.95 


(b) 


9.58-11.34 


10.87-13.39 


13.53-16.82 


26.03-29.23 


14.25-16.75 


64.93-67.16 






0.97-1.98 






39.55-41.50 








2.53-3.00 






17.98-19.12 


(c) 


9.18-10.63 


10.88-12.43 


12.76-15.55 


25.96-29.17 


14.35-17.40 


64.53-67.01 






0.89-1.85 






39.27-41.37 








2.53-2.85 






17.83-19.05 


(d) 


9.57-12.74 


10.84-14.91 


13.50-19.81 


26.18-29.80 


14.19-17.46 


65.05-68.44 



1.02-2.59 39.70-42.44 
2.54-3.08 18.06-19.62 



68%, 39%-42% and 18%-20% for SRCl, SRC2 and 
SRC3, respectively. To understand the behaviour 
of SRC, we plot in Fig. 1 the neutrino poten- 
tial Hn (r, A)=ij|'^^ (r. A) p{r) with three differ- 
ent parameterizations of the SRC. The potential in- 
cluding only FNS is peaked at origin whereas the 
peaks due to F-l-SRCl, F-I-SRC2 and F-HSRC3 are 
at r « 0.8 fm, 0.7 fm and 0.5 fm, respectively. The 
visible reduction of the area under the curves is the 
main cause behind the large changes reported in 
Table III. 




FIG. 1: Radial dependence of Hn (r, A)= (r, A) f{r) 

for the three different parameterizations of the SRC. In the 
case of FNS, /(r) = 1. 

(vi) The maximum variations in M*^"''^ {M^^-^^^) due to 
PQQHHl, PQQ2 and PQQHH2 parameteriza- 
tions with respect to PQQl interaction, but for the 
pathological case "°Ba, are about 21.0% (13.0%), 
19.0% (17.0%) and 18.0% (27.0%). 

(vii) The effect of deformation on A/*^^^ is quantified by 



the quantity D^^) as the ratio of M^^) at zero de- 
formation {(qg — 0) and full deformation {(gg = 1) 
and is given by ^] 



M(^){Cgg^O) 

A/W(C,, = 1) 



(16) 



In Table llVi we tabulate the values of Z?^^' due 
to exchange of light and heavy neutrinos for ^^Ru, 
i02pd, ^°^Cd, i24xe, i30Ba and ^^^Dy nuclei. It 
is observed that, due to deformation effects, the 
NTMEs Af(^) are suppressed by factor of 1.7-10.8 
in the mass range A = 96 — 156. Thus, the de- 
formation plays a crucial role in the iP'^ /3~^)q^ and 
(e/3+)Qj^ modes. 

(viii) It is also observed that excluding the pathological 
case i30Ba the ratios of NTMEs M'^^"^ / A/(°^) are 
about 21-25, 32-37 and 42-49 for SRCl, SRC2 and 
SRC3, respectively. The spread in the above men- 
tioned ratios increases to 21-29, 32-43 and 42-57 
with the consideration of ^^°Ba isotope. 



B. Radial evolution of NTMEs 

In the Majorana neutrino mass mechanism, the radial 
evolution of NTMEs can be studied by defining 



(17) 



The study of radial evolution of NTMEs A/^"") in the 
QRPA by Simkovic et al. "42], and in the ISM by 
Menendez et al. [58|, has established that the magni- 
tude of C'^'^"^ for all nuclei undergoing (/3~/3~)oi, decay 
exhibit a maximum at about the internucleon distance 
r « 1 fm, and that the contributions of decaying pairs 
coupled to J = and J > almost cancel out beyond 
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TABLE 
tion. 



IV: Deformation ratios (i) D'""^' and {ii) D^^^^'> of 
and {eP^)^^ modes for the PQQl parameteriza- 



Nuclei 




F 




F+SRC 






F+SRCl 


F+SRC2 


F+SRC3 






2.53 


2.54 


2.53 


2.53 




(n) 


2.44 


2.34 


2.40 


2.43 


102pj 


(*) 


2.68 


2.74 


2.68 


2.67 




(ii) 


2.40 


2.45 


2.42 


2.40 


lo^Cd 


(i) 


1.92 


1.96 


1.93 


1.92 




(^^) 


1.72 


1.77 


1.74 


1.73 


i^^Xe 


(*) 


3.82 


3.92 


3.84 


3.82 




(n) 


3.36 


3.50 


3.42 


3.38 


i^OBa 


(^) 


4.61 


4.75 


4.65 


4.61 




{^^) 


3.97 


4.15 


4.06 


4.01 




(i) 


10.78 


10.82 


10.76 


10.75 




(^^) 


10.25 


10.01 


10.14 


10.20 



r w 3 fm. In the PHFB model, the radial evolution of 
NTMEs M(°'') and M*^"^^ for (/3~/3")o^ decay due to the 
exchange of light [35] and heavy Majorana neutrinos [IBl 
has also been studied and similar observations have been 
reported. 

Presently, we study the radial dependence of C^^"^' as 
well as C^"^) for i^+(3+)o^ and {eP+)ou modes of ^^^Ru, 
i02pd, loe^d, i24xe, i30Ba and ^^^Dy isotopes in four 
cases, namely F, F+SRCl, F+SRC2 and F+SRC3. In 
Fig. [1 we plot the radial dependence of C^o^^^ and C(°^' 
for ^°^Cd, employing the PQQl parameterization of the 
effective two body interaction, for four combinations of 
FNS and SRC. In Fig. 3, the radial evolution of C'^^"^ 
and C^^^'^ are displayed together for the six nuclei under 
study, for the four combinations of FNS and SRC. 

In the case of light Majorana neutrino exchange, it is 
noticed that the C^"'^-' are peaked at r = 1.0 fm for finite 
size nucleons and the addition of SRCl and SRC2 shifts 
the peak to 1.25 fm. However, the position of the peak 
remains unchanged at r — 1.0 fm with the inclusion of 
SRC3. The radial distributions of C^"'') extends up to 10 
fm although the maximum contribution to M^"") results 
from the distribution up to 3 fm. In the case of heavy 
Majorana neutrino exchange, the C^^^'' are peaked at 
r « 0.5 fm in the case of FNS, and with the addition of 
SRCl and SRC2, the peak shifts to about 0.8 fm, and to 
0.7 fm for SRC3. The radial distributions of C(°^) extend 
up to 2 fm and the total distribution contributes to the 
evolution of M^^^\ Remarkably, the above observations 
also remain valid with the other three parameterizations 
of the effective two-body interaction. 



C. Uncertainties in nuclear transition matrix 
elements and nuclear sensitivity 

The uncertainties associated with the NTMEs M'^^"^ 
and M(o^) for (/3+^+)oz. and (e^+)o,. modes of ^^Ru, 
i02pd, ^°^Cd, i^-^Xe, "OBa and ^^^Dy isotopes due to 
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FIG. 2: Radial dependence of C^°''^(r) and C^°^'(r) for the 
and (e/3+)o^ modes of "''Cd isotope. 



the exchange of light and heavy neutrinos, respectively 
are evaluated by calculating the mean and standard de- 
viation given by 



M 



[K) 



TV 



[K) 



and 



AM 



1 



N 



M 



(K) 



M 



1/2 



(18) 



(19) 



The twelve NTMEs due to the exchange of light as well 
as heavy Majorana neutrinos listed in the three columns 
4-6 and 11-13 (F-l-S) of TableUare employed in this sta- 
tistical analysis for the bare and quenched values of axial 
vector coupling constant gA — 1.254 and gA = 1-0, re- 
spectively. Further, the effect due to the Miller-Spenser 
parameterization of Jastrow type of SRC is estimated by 

evaluating the same mean m''^^ and their standard devi- 
ations AM^^^ for eight NTMEs calculated using SRC2 
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FIG. 3: Radial dependence of C^°''\r) (left )and C(°^'(r) (right) for the (/3+/3+)q^ and (£/3+)o^ modes of ^®Ru, ^^^Pd, 
i"«Cd, i^^Xc, i^'Ofia and ^^®Dy isotopes. In this Fig., (a), (b), (c) and (d) correspond to F, F+SRCl" F+SRC2 and F+SRC3, 
respectively. 
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TABLE V: Average NTMEs m'^' and uncertainties AM^''"'^ for and (e/3+)o^ modes of ^'^Ru, "^Pd, ^"''Cd, ^^''Xe, 

^^"Ba and ^^^Dy isotopes. Both bare and quenched values of qa axe considered. Case I and Case II denote calculations with 
and without SRCl, respectively. In column 9, (1) and (s) denote large and small basis, respectively. 



Nuclei QA Light neutrino exchange Heavy neutrino exchange 

Case! Case II QRPA QRPA SQRPA MCM Case! Case II 















[23] 


[59] 


[60] 






AI?^"^^' 






*Ru 


1.254 


4.59 


0.34 


4.82 


0.11 


3.60 


4.228 




2.383 


148.13 


50.27 


178.29 


29.19 




1.0 


5.13 


0.40 


5.39 


0.13 










165.91 


59.74 


201.49 


35.61 


102pj 


1.254 


4.71 


0.60 


4.97 


0.50 










160.83 


58.13 


195.01 


35.86 




1.0 


5.34 


0.71 


5.63 


0.59 










181.01 


69.02 


221.36 


43.39 


cd 


1.254 


7.57 


0.89 


7.97 


0.72 


4.56 


4.778 


7.85(1) 


3.394 


249.89 


89.73 


302.93 


54.54 




1.0 


8.52 


1.04 


8.98 


0.84 






8.17(s) 




281.22 


106.57 


343.82 


66.13 


i24Xe 


1.254 


3.50 


0.42 


3.69 


0.32 


5.27 


2.975 




8.301 


124.84 


44.75 


151.45 


26.71 




1.0 


3.96 


0.49 


4.19 


0.37 










140.70 


53.16 


172.11 


32.44 


i^Ofia 


1.254 


2.60 


0.80 


2.75 


0.82 


5.52 


5.579 




5.130 


97.35 


41.65 


118.11 


34.03 




1.0 


2.94 


0.90 


3.12 


0.92 










109.75 


48.71 


134.24 


39.56 


i5«Dy 


1.254 


2.22 


0.31 


2.33 


0.29 










72.96 


26.36 


87.78 


18.02 




1.0 


2.50 


0.36 


2.63 


0.33 










81.66 


31.12 


99.13 


21.41 



TABLE VI: Limits on the effective mass of light < > and heavy < Mjv > Majorana neutrinos for the (/3^/3^)q^ and 
(e/3+)jj^ modes of ''^Ru, ^^^pd, lOScd, ^^^Xe, ^'^°Ba and ^^'^Dy isotopes. 





1/2 


(yr) 


Ref. 




(m 


.) (eV) 


(A/iv> 


(CeV) 


emitters 










/?+/?+ 




/?+/?+ 




^•^Ru 


>3.1xl0i'' 


>6.7xl0''' 


[62j 


1.254 
1.0 


4.02 xlO'' 
5.66x10^ 


7.94 xlO'' 
1.12x10^^ 


4.41x10 
3.17x10 


2.23x10^ 
1.61x10^ 


i"«Cd 


>1. 2x10^1 


>2. 2x10^1 


m 


1.254 
1.0 


1.16x10^ 
1.62x10^ 


2.27x10^ 
3.16x10^ 


1.57x10*1 
1.13x10*1 


8.04x101 
5.80xl0i 


i24xe 


>4.2xl0i'' 


>1.2xl0i^ 


[63] 


1.254 
1.0 


1.22x10^ 
1.70x10^ 


1.68x10'' 
2.33x10*1 


1.61x10^ 
1.16x10^ 


1.17x10^ 
8.46x10^ 


""Ba 


>4.0xl0^i 


>4.0xl0^i 


m 


1.254 
1.0 


4.11x10^ 
5.70x10^ 


4.19x10^ 
5.82x10^ 


5.01x10^ 
3.62x10^ 


4.91xl0i 
3.55xl0i 



and SRC3 parameterizations. In Table |Vl we display 
the calculated averages and their variances along with 
all the available theoretical results in other models for 
96Ru, i02pd, i06cd^ i24xe, i30Ba and ^^^Dy isotopes. 

In the case of light Majorana neutrino exchange, it is 

observed that the uncertainties AM^^'^^ but for ^^°Ba are 
about 7%-14% and the exclusion of NTMEs M'^^"^ cal- 
culated with the Miller-Spencer parameterization of Jas- 
trow SRC, reduces the uncertainties to 2%-12% for both 
QA = 1.254 and gA — 1-0. Pathologically, the uncertainty 

AM' ' « 30% in the case of ^"^"Ba, remain unaltered 
due to the large effects of PQQHH2 parameterization. 

The estimated uncertainties AM^^'^'' but for ^•^°Ba iso- 
tope in the heavy Majorana neutrino mass mechanism, 
are about 35% for gA = 1.254 and gA = 1.0. Estimation 



of uncertainties for eight NTMEs M calculated using 
the SRC2 and SRC3 parameterizations again reveal that 

the AM*""^** are reduced to 16%-21% due to the exclu- 
sion of SRCl. In ^-^^Ba isotope, the same pathological 
behaviour is noticed. 

In the QRPA calculations of Hirsch et al. [l^j and 
Staudt et al. [Hi], the NTMEs M^^"^ are almost identi- 
cal but for "'^^^Xe, in which the difference is approximately 
by a factor of 1.8. Stoica et al. [60] have used SQRPA 
model with two model spaces, namely small basis (oscil- 
lator shells of ifiuj — bhuj -\- 113/2 orbit) and a large basis 
(oscillator shells of 2'huj — 5fiu! + ii3/2 orbit). They used 
the same SPEs as those of Hirsch et al. and an effective 
two-body interaction derived from the Bonn-A potential. 
The NTMEs calculated in the SQRPA W\ do not depend 
much on the model space and differ by a factor of 1.8 ap- 
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TABLE VII: Predicted half-lives, corresponding extracted effective mass of heavy Majorana neutrino < A/jv >, nuclear sensi- 
tivities and ^'"^^ due to exchange of light and heavy neutrino, respectively, for {P^P^)^ and (e/3+)j| modes of ^^Ru, 
i02pd, lo^cd, i24Xe, i30Ba and ^^"Dy isotopes. 



e+pp 






= 0.05 eV) 


{Mm) (GeV) 








(ON) 


emitters 






eP+ 




P+P+ 


eP+ 


P+P+ 


eP+ 


"'^Ru 


1.254 
1.0 


2.01+H^xlO^o 
3.97l°;i^xl03o 


1.69+0'0«xl0^" 


3.55+JJ;^»xlO» 
3.59lo;e3XlO* 


0.721 
0.513 


2.486 
1.768 


26.70 
19.19 


92.02 
66.13 


102p^ 


1.254 
1.0 


- 


7.04+i|^x 10=^0 


3.76+r6'^xl0«t 
3.77lO;74xio8t 


- 


0.385 
0.277 


- 


15.11 
10.90 


Cd 


1.254 
1.0 


6.4911:3^ X1Q29 
1.27lE!;2iXl0^° 


4.52l0;?txl02« 


3.64l0;f^xl0« 
3.67l!5;7ixlO* 


1.268 
0.908 


4.806 
3.442 


48.20 
34.79 


182.6 
131.8 


i24xe 


1.254 
1.0 


2.51+J;-^^x1Q30 
4.841100x103° 


1.35+0^1 xlO^^ 


3.93+^69 xlO« 
3.94lo;^^10« 


0.645 
0.465 


2.777 
2.001 


26.44 
19.11 


113.9 
82.30 


"ORa 


1.254 
1.0 


2.70j:?'^^x 10=^1 
5.20l^;?ix 10=^1 


2.82+^*1 X 10^^ 
5.42l^j^xl02« 


4.1211;!^ xlO® 
4.13li;22xl0* 


0.197 
0.142 


1.926 
1.387 


8.45 
6.11 


82.68 
59.76 




1.254 
1.0 




5.94+i'^;^xl029 
1.1510^1x10^0 


3.62+0'^txl0«t 
3.6llo>8xlO«t 




1.327 
0.953 




50.04 
35.94 



"t denotes (£/3+)o.y mode only. 



proximately from those of Hirsch et al. [23| . 

There are no available theoretical results and experi- 
mental half-life limits for the ^°^Pd and ^^^Dy isotopes. 
The extracted limits on the effective light neutrino mass 

< nil, > as well as heavy neutrino mass < > using 
the phase space factors given in Ref. [13] and presently 
available experimental limits on observed half-lives of 
{(3'^ (3'^)qi, and {£(3'^)qj, modes are presented in Table IVll 
The extracted limits on < rrii, > and < Mn > are 
not so much stringent as in the case of (/3~/3~)oj^ de- 
cay. Moreover, better limits are obtained in the case 
of (e/3+)Qj^ mode even for equal limits on half-lives of 
{P'^P~^)qi, and (e/3''')oi/ niodes. The best obtained limits 
for lo^cd isotope are < > < 1.16 x 10^ (2.27 x 10^) 
eV and < Mn » 1.57 x 10* (8.04 x 10^) GeV in case of 
(/3"'"/3'*')q^ and (e/3^)o,y modes, respectively. In the case 
of (/3^/3^)o^ and (£/?'*')oy modes, the extracted limits on 
the effective neutrino masses < > and < Mjv > 
are not stringent enough and hence, we calculate half- 
lives of these modes to be useful in the design of fu- 
ture experimental setups. The half-lives of (/?"'' /3^)o^ and 
(£/3+)qj^ modes for < nii, >= 50 meV are calculated and 
extracted corresponding limits on heavy neutrino mass, 

< Mjv >, are given in the same Table IVlIl 

In the absence of stringent limits on the effective neu- 
trino masses < my > and < Mn >, it is useful to calcu- 
late the nuclear sensitivity, defined as [1^ 



(20) 



where K stands for Oi> or OA'^ mode and an arbitrary 
normalization factor 10^ is introduced so that the nuclear 
sensitivity turns out to be order of unity. 



It is observed that in general, nuclear sensitivities 
for (£/3^)qj, mode are larger than those of (/3^/3^)Qy 
mode. Further, the nuclear sensitivities for (/3^/3^)Qy 
and (e/3+)o, modes of i06cd, ^^Ru (^^^Xe), ^^^Xe (^^Ru), 
^^°Ba, ^^^Dy and ^^^Pd isotopes, respectively, are in the 
decreasing order of their magnitudes. 



IV. CONCLUSIONS 

We have calculated sets of twelve NTMEs M^""^ and 
M(o^) for (/3+/3+)o^ and (e/3+)oy modes of ^^Ru, ^^^Pd, 
i°^Cd, i^^Xe, ^^°Ba and ^^^Dy isotopes by employing 
the PHFB model with four different parameterizations 
of the pairing plus multipolar type of effective two body 
interaction and three different parameterizations of the 
short range correlations. To estimate statistically the 
uncertainties in NTMEs, mean and standard deviations 
of sets of twelve NTMEs M^"^) and M(°^) calculated 
with dipole form factor and short range correlations are 
employed for both gA = 1.254 and gA = 1-0. It is ob- 
served that the largest standard deviation turns out to 
be around 30% in the case of ^'^'^Ba isotope due to the 
dominant contribution of deformation in PQQHH2 pa- 
rameterization. But for ^■^'^Ba, the maximum uncertainty 
in NTMEs M^^"^ is around 14%, which becomes smaller 
by 2% excluding the NTMEs calculated with SRCl in the 
case of ^^®Dy isotope. The uncertainties in M'"^^ due to 
the exchange of heavy Majorana neutrino is about 35%. 
Exclusion of NTMEs calculated with SRCl, reduced the 
uncertainties by 14%-19%. 
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